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We numerically investigate a bosonic representation for 
hole pairs on a two-leg t-J ladder where hard core bosons on 
a chain represent the hole pairs on the ladder. The interaction 
between hole pairs is obtained by fitting the density profile ob- 
tained with the effective model to the one obtained with the 
t-J model, taking into account the inner structure of the hole 
pair given by the hole-hole correlation function. For these in- 
teractions we calculate the Luttinger liquid parameter, which 
takes the universal value K p = f as half filling is approached, 
for values of the rung exchange J' between strong coupling 
and the isotropic case. The long distance behavior of the hole- 
hole correlation function is also investigated. Starting from 
large J', the correlation length first increases as expected, but 
diminishes significantly as J' is reduced and bound holes sit 
mainly on adjacent rungs. As the isotropic case is approached, 
the correlation length increases again. This effect is related 
to the different kind of bonds in the region between the two 
holes of a hole pair when they move apart. 



I. INTRODUCTION 

The study of strongly correlated electrons in ladder 
systems has been an active field in recent years. Ladders 
consisting of a number of coupled legs (or chains) show a 
rich variety of phases (l]-(|] depending on the number of 
legs and the electron density. One property of particular 
interest is the binding of holes into pairs when weakly 
doped into a half-filled two leg ladder, to form a Luther- 
Emery liquid. In this liquid the spin degrees of freedom 
are gapped and the low energy sector can be mapped to 
a boson liquid of hole pairs. This can be studied an- 
alytically in the weak coupling limit described, for ex- 
ample, by a one-band Hubbard model with weak interac- 
tions. In the strongly interacting limit described by a t-J 
model, 0| one must use numerical methods such as Lanc- 
zos diagonalization and the recently developed Density 
Matrix Renormalization Group (DMRG) method. [|]|] 
Large simulations can be carried out using the DMRG 
yielding information on the equal time correlations. 

In this paper our goal is to use the hole density dis- 
tributions determined by DMRG to determine the form 
of the effective boson model that we know on general 
grounds describes the hole pairs. We use these DMRG 
results to parameterize their dispersion and their effec- 
tive mutual interactions and also the interaction with 
hard wall boundary conditions that occur in DMRG. In 
this first paper we limit ourselves to the simplest case of 
the two leg ladder. The same approach can be made for 



wider ladders with an even number of legs and is cur- 
rently in progress. 

We start with the t-J ladder Hamiltonian 

H = ~t ]T V(4^C l+1 ^ a + 4 fl)J >C i J,<r)^ (1) 

+J' y^(Sj,iS»,2 - jni^m^) 

i 

where i runs over L rungs, j(= 1,2) and a (=T>4) are 
leg and spin indices. The projection operator V = 
Eli j0-~ n id,1 n i,j,l) prohibits double occupancy of a site. 
Unless noted otherwise we set t — t' and J = 0.35 1. 



II. THE STRUCTURE OF A HOLE PAIR 

In this section we examine a hole pair on a t-J ladder, 
studying the internal hole-hole correlation (hhc) func- 
tion. The hhc function gives us information about the 
inner structure of a hole pair and depends on the param- 
eters t, t' , J and J' in the t-J Hamiltonian given above. 
Since the DMRG computations are performed with open 
boundaries, special care must be taken in the measure- 
ment of the hhc-function, which we define as 



9j(ri,r 2 ) 



(2) 



The index j — 1,2 denotes intra- and inter- leg- 
correlations respectively, rtfj = 1 — denotes the hole 
density at site and n h the average hole density. In- 

troducing relative and center of mass coordinates, r and 
i?o, we measure gj(r,R ) on (40 x 2)-ladders with open 
boundaries for Rq ~ L/2. In a range Rq = L/2 ± 5, 
gj(r,Ro) depends only weakly on Rq. Furthermore, we 
have tested that better accuracy for the hhc-function can- 
not be achieved on larger ladders up to (60 x 2) sites. We 
conclude that these measurements give good values for 
gj(r) on an infinite ladder. 
We denote by 



(3) 
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the sum of intra- and inter-leg correlations. We have ex- 
amined g(r) as a function of the rung exchange coupling 
J'. As can be seen in Fig. |], a significant change in the 
hhc-function occurs when the isotropic case (</' — > J) 
is approached from the strong rung coupling regime 
(J' ^S> J). In the latter case, the holes sit mainly on 
the same rung whereas in the former case, they are found 
predominantly on adjacent rungs. We study the long dis- 
tance behavior of the hhc-function by fitting the tail to an 
exponentially decaying function, g(r) ~ exp(— r/^hhc)- 
The results are shown in Fig. |2|, where ^hhc is plotted 
as a function of J' . Starting from strong rung coupling, 
the correlation length first increases as expected, but de- 
creases significantly as J' is reduced and bound holes sit 
mainly on adjacent rungs. As the isotropic case is ap- 
proached, the correlation length increases again. 

This behavior can be explained by looking at the dif- 
ferent kinds of bonds in the region between two holes of 
a pair as the holes move apart. In Fig. ||, we show the ex- 
change fields around a pair of holes. These are obtained 
by measuring Sj ■ Sj after projecting out a particular con- 
figuration of two holes from the ground state \ip) given 
for two holes on a (16 x 2)-ladder. Denoting the cor- 
responding projection operator by Ph, Ph Si • Sj Ph has 
been measured and normalized by (ip\Ph\ip). The proce- 
dure is well described by White and Scalapino. [|10| This 
kind of measurement gives us a "snapshot" of the spin 
configuration around a dynamic hole. In the following wc 
use the term "bond" simply to indicate that ( Sj-Sj ) < 0. 
If this expectation value is close to —0.75 for two sites i 
and j, we say that there is a "singlet bond" connecting i 
and j. 

For J' > 1.2 1 the rung exchange dominates. When a 
hole pair virtually splits up, the region between the two 
holes reverts to rung singlets, as shown in Fig. || A high 
energy is needed to split up the holes, but after they are 
separated there is no further increase in energy in moving 
them apart, since the diagonal exchange bonds are weak. 
For J' < 0.7 1 the region between the holes does some- 
thing more complicated, which facilitates hole hopping at 
the expense of higher exchange energy on the interven- 
ing rungs. As shown in Fig. |3|, various diagonal exchange 
bonds are formed, which turn into the rung singlets when 
the holes hop back together. This allows for easy hopping 
and low kinetic energy, but the exchange energy cost is 
high. For J' = 0.7 1, more and more distorted exchange 
bonds are created as the holes move apart. This leads to 
a stronger interaction for larger separation than in the 
case for J' = 1.2 1. This weakening of the long range 
interaction between the two holes explains the increase 
of Chhc in Fig- § in the range between J' = 0.7 1 and 
J'~1.2t. 

Somewhat similar analysis can be found in Ref. |l(]] 
where the authors have argued that the exchange bonds 
seen in Fig. ^ are responsible for the pairing seen in the 
isotropic case, as well as for stripe formation. 



III. THE EFFECTIVE MODEL FOR HOLE PAIRS 

In the strong coupling limit J' 3> J, in the ground 
state, two bound holes are on the same rung and a de- 
scription of tightly bound hole pairs moving in a back- 
ground of singlet rungs is appropriate. Considering these 
hole pairs as hard core bosons (heb) one can map to an 
effective boson model as discussed in Ref. Q and obtain 
the hole density directly from the corresponding heb den- 
sity as well as the correlations between hole pairs. 

Here we propose an effective model that applies for 
any value of J' between strong coupling and the isotropic 
case. Since holes on a ladder can pair with more weight 
on adjacent rungs, our effective model should incorpo- 
rate the possibility that the "center of mass" of a hole 
pair can lie on a rung or between two rungs as shown 
in Fig. ||. Note that for even (odd) distance r along the 
legs between two holes, the center of mass lies on a rung 
(between two rungs). 

To motivate our effective model we study first the case 
of pairing between two holes on a two-leg ladder with 
length L where the Hamiltonian is given by 

2 L 

H = E + b Uuhj) (4) 

3=1 i=l 
L 

-t h J2( b li b ^ + b i,2Ki) + Vh ■ 

i=l 

Here b\ and bi denote hole creation and annihilation op- 
erators respectively. Only single occupation of a site is 
allowed and we use periodic boundary conditions. We 
choose the attractive interaction Vh between the two 
holes strong enough to give an even parity bound state. 
We want to find an effective Hamiltonian for the paired 
holes, expressed in terms of a new effective boson oper- 
ator Bj which describes the pair as a whole, and then 
answer the question, how do we obtain the hole density 
n\j = b\ jbij from the density Ni = P\Bi of the effective 
bosons? 

For the Hamiltonian above, with the help of the oper- 
ators 

PR,r = ^= ( b R-Z,l b R+$A + 4-5,2^+5,2) C 1 ~ <M (5) 



PR,r 



V2 

"r-5,i"r+5,2 



r = 



(G) 



the wavefunction for the ground state can be writ- 

ten as 

* t = 7f E E (Mr)pi,r + Mr)pR, r ) (7) 



R— int. r eve n 



TTjf? E E {^( r )PR,r + < P 2 ( r )PR,r) 



fi=ha!f r odd 
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Here R and r(= 0, 1, 2, . . .) denote the center of mass and 
relative coordinates along the legs as sketched in Fig. 
<Pj (r) depends only on r and the interaction Vh in Eq. 
The pair correlation function gj (r) for the holes b\ in this 
model is simply given by gj(r) = \<Pj(r)\ , where j = 1,2 
denote intra- and inter-leg correlations respectively. 

The probability of finding the center of mass of a pair 
centered on a rung, w int , or between two rungs, Whaif , is 
given by 



Mint 



£Mr)| 2 + |^(r)| 2 = $>(r) (8) 



w ha]S = \<Pi Wl 2 + W2(r)\ 2 = ^,9(r) 

Todd ''odd 

where g(r) denotes the sum of the inter- and intra-leg 
correlations as defined in Eq. (||). Note that the proba- 
bility to find the pair on or between rungs depends only 
on whether R is integer or half integer. 

The same occupation probabilities can be obtained 
with the Hamiltonian 

l L+h 



for one boson B R which moves on a closed chain of length 
L under the action of a periodically varying onsite poten- 
tial e. Here N R = B R B R . 

Figure ^ shows the mapping of hole pairs from the 
ladder to effective bosons on a single chain. Note that 
the center of mass coordinate R of the pairs determines 
the position of the effective boson. Here the ratio of the 
probabilities Whaif/ifmt to find the boson on a site with 
half integer or integer R depends only on e/t* and can 
easily be obtained as 



w hal{ _ e 2 + (4i*) 2 - e^+Wf 



Mint 



+ {At*) 2 + ey/e 2 + (4i*) 2 



(10) 



From this equation and Eq. (ph we obtain for e 



e = 4 1* sinh (In 



Er„d d 9(r) 



(11) 



The boson operator B R can be expressed in terms of 
the operators p' R r and p R r as 
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(<PlO) PR,r + f2{r)p R<r ) 



g( r ) 



for int. R 



for half int. R. 



(12) 



Once the density Nr for the model (||) has been com- 



puted, we obtain the density n 



T-,3 



f° r the model 



([|), taking into account the inner structure of the effec- 
tive boson (n2h by the convolution 



1 Er_ 



N„ 



i+i 



, lE rodd «7(r)(JVU 5 +iVi + i) 



(13) 



We are interested in the ground state properties of 
the t-J model. Since the spin part of the ground state 
wavefunction is antisymmetric and the spin excitations 
are gapped, the charge degrees of freedom can be 
described considering hole pairs as effective (hard core) 
bosons moving in a spin liquid. 

To obtain the effective Hamiltonian describing hole 
pairs in the t-J model (0), we model the holes as in 
Eq. (Eh and define in analogy to the preceding 



L—i 



H = -t*J2 (B R B R+ i+B R+h B R ) 
+e J2 N R + y int + V b 



(14) 



where N R = B R B R and the B R and B R denote heb 
creation and destruction operators respectively. Since 
we are using open boundaries, we have to take into ac- 
count the interaction of the hole-pairs with the bound- 
aries. The potential Vb has been introduced to describe 
this effect. The potential V5 n t gives the interaction be- 
tween heb, i.e. hole pairs in the t-J model. The onsite 
potential e on half-integer sites has the same meaning as 
explained above. 



IV. COMPUTING THE MODEL PARAMETERS 

In this section we determine the effective model pa- 
rameters t* , e and the interactions Vb and Vint- 



A. Onsite potential e 

The onsite potential e has been calculated with the 
hhc- function g(r) obtained numerically by the DMRG 
for two holes on (40 x 2) t-J ladders for various J'. The 
results are given in Fig. They clearly show that for the 
isotropic case, i.e. J = J' , where e < 0, the hole pair is 
mainly centered between two rungs, whereas for strong 
coupling, i.e. J' S> J, where e > 0, both holes of a pair 
sit on the same rung and can be well described by the 
effective model given in Ref. ||] . 
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B. Interaction with the boundaries 

As e is determined by the hhc-function alone, we can 
compute the hcb density JV», convolute it with g(r) ac- 
cording to Eq. ( p^| ) and compare it with the hole density 
of the corresponding t-J system. The hole density n^j 
can be obtained using Eq. ( |l3| ) by identifying n\j with 
n^j. In this way we obtain Vb and Vj n t by fitting the 
density profile obtained with the effective model to the 
one obtained with the t-J model. Since Vint in Eq. (^) 
gives no contribution for one hcb, we can obtain Vb by 
considering one hole pair in the t-J model. We choose 
an exponentially decreasing form for Vb 



V b = v b N R (« 



(15) 



We have tested the results obtained in this way by 
comparing the density profiles for various numbers of hole 
pairs and system lengths. From strong coupling down 
to the isotropic case we find good agreement between 
the density profiles obtained from the t-J model and the 
effective model, as can be seen in Fig. [[o| 

We should note here that the fits for Dj nt were not 
as stable as the fits for lib- Varying simultaneously the 
parameters v\ , v and £ in Eq. ( |l5| ) , one can find another 
set of parameters which also gives a good fit. However, 
we found that the Luttinger liquid parameter K p , which 
we calculate in the next section, only weakly depends on 
this ambiguity. Using only two parameters to determine 
Dint suppressing v\ in Eq. (|l5|), one can also fit the density 
profiles with a different set of parameters, but in this case 
K p is strongly affected. 



with increasing distance from the boundary. Figure ^ 
shows fits for J' — 0.35 1 and J' = 10 1. and the param- 
eters Vb and are shown in Fig. |?]. Up to J' — 3.0 1, Vb 
is positive and £b is finite. Above this value u& becomes 
negative and £b is zero. In this case only the outermost 
rungs become attractive for hole pairs as can be expected 
for large J', due to the charge part 
J' term. 



3 n i,j n i+hj of the 



C. The interaction term Vint 

To obtain the interaction potential Vi nt we proceed in 
the same way as for Vb- We choose a hard core analytical 
form 

v ** =EE v ^(\ R - N R N R' ( 16 ) 

R R'>R 

oo if r < r min 

v e « ~ if r > r mm 

with r m i n > 1 and considered (30 x 2) t-J ladders with 
four holes and the corresponding effective model with 
two hcb. The infinite repulsion for r = \ reflects the 
hard core condition for hole pairs, so that the minimum 
distance of two hole pairs in the t-J model is 1. 

Since the hole pair becomes broader with decreasing 
J', it was useful to let r m i n vary. Again, we have made 
fits between the density profiles as shown in Fig. || and 
obtained the parameters represented in Table |. Fig- 
ure ^| shows Di n t for different J 1 values. Starting from the 
isotropic case, the interaction is repulsive and becomes 
less long ranged as J' increases. For J' = 5i, v- mt be- 
comes negative for r = 2. The same holds for J' = 10 1, 
where v- mt vanishes for r > 2. This can be expected from 
a perturbative approach as can be found in Ref . [pj , where 
the hole pairs on a t-J ladder have been mapped to hcb 
on a chain with next-nearest neighbor interaction. 



D. The hopping matrix element t* 

Up to now the hopping matrix element t* did not play 
any role. We were only interested in the hcb density 
profile, which is not affected by the energy scale fixed by 
t*. 

We compute t* by finite size scaling with the ansatz 
that the difference between the ground state energy for 
two holes i?2h and for zero holes Eq^ in the t-J model is 
given by 



E 2h (L) - E oh (L) ~ const + t* E oS {L) 



(17) 



Here E e g denotes the ground state energy of the corre- 
sponding effective model with one boson and with t* = 1. 
We obtained t* by fitting Eq. ( p7| ) to a straight line. The 
result is shown in Fig. |ll|. 

According to the effective model given in Ref. ||, one 
can obtain the corresponding effective hopping matrix 
element t° by a finite size scaling with the ansatz 



Ezn(L) - E oh (L) ~ const + <°(tt/(L + l)) 5 



(18) 



Here, E2h(L) and Eo^(L) are for the t-J model. A per- 
turbative estimate for t° is given in Ref. || . For large J' 
one obtains 



t° = 2t 2 /(J' -4t 2 /J') 



(19) 



Figure p"T| shows t° versus J'. It can be seen that for J' — > 

00 the t° obtained by perturbation theory tends to the 
one obtained by finite size scaling. Note that t° stands 
for the hopping of a hole pair from one rung to another, 
whereas in our effective model t* mediates between i and 

1 + | . We can obtain t° from t* by expanding the single 
boson dispersion of the effective model for small momenta 
k. In this way we obtain 



t° 



t* 



^( 2 + {At*f 



(20) 



4 



Figure |ll| shows that t° obtained in this way from our ef- 
fective model fits very well to the one obtained by finite 
size scaling. 

An independent test for t* comes from the comparison 
of the inverse compressibility kT 1 obtained from the t-J 
model and the effective model. The inverse compressibil- 
ity could be obtained from 



_ x p 2 d 2 E Q (p) 

1 ip)= n^p^- 



(21) 



Here the volume f2 = 2 L for the t-J model and £1 — 
2 L — 1 for the effective model. With N we denote the 
hole or hcb number and with p = N/Q the corresponding 
density. The second derivative of the ground state energy 
with respect to the particle number N was calculated by 
the discretization 



d 2 E E {N + A) - 2 E {N) + E (N - A) 



(22) 



ON 2 ~ A 2 
Here A = 2 and A = 1 has been used for the t-J model 



and the effective model respectively. Figure 12 shows the 
results for k _1 computed for a (40 x 2) t-J ladder with 2 
holes and 4 holes. As can be expected, the results agree 
better for large J', where the hole pairs are narrower and 
the interaction between pairs is less long-ranged, but the 
agreement is satisfactory for all values of J'/t. 



V. LUTTINGER LIQUID PARAMETER K p 

In the previous section we obtained the interaction be- 
tween two hcb in the effective model (i.e. two hole pairs 
in the t-J model). In order to obtain information about 
the long range correlations we have calculated the Lut- 
tinger liquid parameter K p in the effective model. 

Our effective model belongs to the same universality 
class as the hcb model with next-nearest neighbor inter- 
action, which in turn is equivalent to the XXZ model 
and which has been solved exactly by a bosonization ap- 
proach and conformal field theory. 11 1^] From Ref. Q 
we obtain a power law decay at large distances for the 
charge density wave correlations and the superconduct- 
ing correlations given as 



(N r N ) ~ const x r~ 2 + const x cos(27rpr) r~ 2 Kp (23) 

(BlB ) ~ const x . (24) 

Here p denotes the electron density on the ladder and half 
filling corresponds to p = 1. These relations show that 
the superconducting correlations (BJ,Bq) are dominant if 
K p > i. 

For hcb in one dimension, K p can also be obtained 
from the relations 11131] 



K p = nv c L 



v r K„ 



d 2 E 
ON 2 
7T d 2 E Q (<S>) 
L d<f> 2 



(25) 



Here Eq denotes the ground state energy for a closed ring 
of length L with Nf, hcb and Eq(^>) is the ground state en- 
ergy of the system penetrated by a magnetic flux <f> which 
modifies the hopping by the usual Peierls phase factor, 
1 i— ► £ cxp(± i$>/L). From these two equations the charge 
velocity v c can be eliminated. The second derivative of 
the ground state energy with respect to the particle num- 
ber was calculated by the discretization given in Eq. ( g^ ) 
with A = 1. We used exact diagonalization for system 
lengths between 32 and 220 and with iVj, between 2 and 
4. Analogously, the second derivative with respect to the 
magnetic flux was calculated for the same configurations. 

We calculated the Luttinger liquid parameter K p for 
the interaction potentials given by the parameters in Ta- 
ble | and found the universal value K p = 1 as half fill- 
ing is approached, for any value of J' between strong 
coupling and the isotropic case (Figs. |l3| and [l4|). For 
Nb/L — > 0, corresponding to a very dilute hcb gas, we 
have K p = 1 + 0(Nb/L), independent of the value of 
the interaction and consistent with Refs. |I^,|l5|. Down 
to p ~ 0.875 the superconducting correlations are domi- 
nant, since K„ > \ and since we are far away from phase 
separation. || Below p ~ 0.875, K p becomes less than 
one half and at one quarter doping (p — 0.75) K p < | 
for J' < t. Figure [13 shows also earlier results on small 
clusters from Ref. |^ and Ref. [19], obtained using ex- 
act diagonalization. The deviations from our results are 
most probably due to finite size effects close to half filling 
in Ref. § and Ref. (l6). 

We briefly want to discuss the possibility of commen- 
surability effects. At a filling of p = 0.75, one might ex- 
pect that commensurability effects would stabilize long 
ranged charge density wave ordering, since K p is quite 
small (K p — 0.232 for the isotropic case). In the finite, 
open systems studied here, a static CDW is pinned by the 
boundaries, as can be seen in Fig. [if] for the t-J model as 
well as in Fig. [IB] for the effective model. A careful analy- 
sis using large systems and finite size scaling is necessary 
to determine if the CDW order is long ranged. The decay 
of CDW oscillations away from the edges of the system 
is quite slow and is consistent with long ranged CDW 
order, but a very slow algebraic decay cannot be ruled 
out. A more detailed analysis of the effective model (us- 
ing 200, 240, 280, 320 and 400 sites) shows a decreasing 
charge gap with increasing system size. Finite size scaling 
is consistent with a vanishing or possibly a small charge 
gap in the infinite system. 



VI. CONCLUSIONS 

The effective model derived in this paper works well 
for the hole density of a two-leg t-J ladder for various 
fillings. In other words, the low energy physics of hole 
pairs on a ladder can be well described by a model of 
hard core bosons on a chain with each boson represent- 
ing a pair of holes. The interaction between the hard core 
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bosons was determined by fitting the density profile ob- 
tained with the effective model to that of the t-J model, 
taking into account the inner structure of the hole pair 
given by the hole-hole correlation function. Starting from 
the isotropic case, with equal exchange couplings on the 
rungs and legs the interaction between two hole pairs is 
long ranged and repulsive but becomes attractive and of 
nearest neighbor type when the strong coupling regime is 
approached. The same holds for the interaction of a hole 
pair with the boundaries. We choose a simple form for 
the interaction between the bosons in order to use only a 
few parameters. The results obtained from the effective 
model are insensitive to the specific ansatz used for this 
interaction. The Luttinger liquid parameter K p has been 
calculated for electron densities from p = 0.982 down to 
p = 0.75 (half filling corresponds to p = 1). Down to 
p ~ 0.875 the superconducting correlations are found to 
be dominant and K p — > 1 for p — > 1. For commensu- 
rate filling, p = 0.75, there might be true charge density 
wave ordering and a charge gap. Further investigations 
are necessary to clarify this question. The hopping ma- 
trix element for the bosons in the effective model, which 
allows one to calculate the inverse compressibility, could 
also be determined. Comparing to the t-J model we find 
good agreement for strong rung couplings and only small 
deviation (less than ±10%) near the isotropic case. 

An interesting feature appeared in the hole-hole cor- 
relation function. The correlation length, £hhc, does not 
monotonically increase as one approaches the isotropic 
case from strong rung coupling. Instead, it decreases in 
the interval 0.7 t < J' < 1.2 t. This unexpected behavior 
can be traced to the fact that the interaction between 
holes is dominated by the simple rung exchange bonds in 
the strong rung coupling regime and by the diagonal rung 
exchange bonds near the isotropic case, J' — > J. Reflec- 
tions of the structure in ^hhc as a function of J', shown 
in Fig. |^, can also be found in other properties, e.g. the 
interaction of the effective bosons with the boundaries 
and the compressibility. 
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FIG. 3. Exchange bonds { S, • Sj ) around two dynamic 
holes after projection: t = t', J = 0.35 1, J' = 1.2 t for the 
upper and J = 0.35 1, J' = 0.70 t for the lower figure. The 
thickness of the lines is proportional to the bond strength. 
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FIG. 1. Hole-hole-correlation function g(r) obtained from 
DMRG calculations on a (40 x 2) t-J ladder for J' = 0.35 1 
and J' — 5.0 t with one hole pair (J = 0.35 1, t — t'). 



TABLE I. Parameters for the interaction potential Vint 
given m Eq. © obtained from the density profiles of a 
(30 x 2) t-J ladder for various J' (J = 0.35 1, t = t'). 



r = odd 



J' 


Tmin 


VI 


V 




0.35 


2 


4.973 


0.478 


0.376 


0.5 


3 
_> 


1.28 


0.580 


0.373 


1. 


1 


10.99 


0.586 


0.439 


2. 


1 


5.00 


0.089 


0.363 


5. 


1 


-0.048 


0.0029 


0.400 


10. 


1 


-0.13 


0.0 






i=l 2 | 

- R = integer t R = half integer 

R = 1 3/2 2 

O O rio ^ e • center of mass pair hcb 

FIG. 4. Mapping of pairs of holes from a ladder to hcb on 
a chain depending on the different positions of their center of 
mass and their relative separation r along the legs. R denotes 
the center of mass coordinate along the legs. 
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J'/t 

FIG. 2. Correlation length ^hhc of the hhc-function g(r) as 
a function of J' (J = 0.35 1, t — t') measured on a (40 x 2) 
t-J ladder. 




FIG. 5. The onsite potential e obtained from Eq. (|ll|) as 
a function of J' (J = 0.35 1, t = t'). The hhc-functions were 
obtained from DMRG calculations on (40 x 2) t-J ladders. 
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0.06 




FIG. 6. Hole density for two holes on a (40 x 2) t-J 
ladder with J' = 0.35 1 and J' = 10* ( J = 0.35 i, t = *') com- 
puted directly and with the effective model. For the effective 
model the data for Vb — and the best fit is shown. 
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FIG. 8. Hole density n t h for J' = 0.35 1 and J' 



calculated for four holes on a (30 x 2) t-J ladder computed 
directly and with the effective model (J = 0.35 1, t = t'). For 
the effective model the data for Vint = and the best fit is 
shown. 
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FIG. 9. Interaction parameter Vi nt (r > r m in) between 
the hcb obtained from (30 x 2) t-J ladders for various J' 
(J = 0.35 1, t = t'). 
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FIG. 11. Hopping matrix elements for our effective model 
and for the model in Ref. || obtained for the large J' limit 
(J = 0.35 t — t'). Explanations are in the text. 




FIG. 10. Hole density n\ for J' = 0.35 t and J" = t for six 
holes on (40 x 2) and four holes on (20 x 2) sites respectively, 
computedwith the interaction potentials obtained from fits 
with two hole-pairs on 30 sites (J = 0.35 1, t — t'). For the 
effective model the data for Vint = is also shown. 




FIG. 12. Inverse compressibility k~ 1 for the t-J Ladder 
computed directly and with the corresponding effective model 
for 2 holes (magnified lOx) and 4 holes for a 40 x 2 t-J ladder 
(J = 0.35t, t = *')■ 
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FIG. 13. K p as a function of p for various J' (half filling 
is given by p = 1, J = 0.35 1, t = i'). 
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FIG. 15. Density profile computed with the effec- 
tive model for a (200 x 2) t-J ladder at one quarter dop- 
ing (p — 0.75) with the parameters of the isotropic case 
(J = J' = 0.35 1, t = t'), showing a pinned CDW. 
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FIG. 14. K p as a function of J' for various electron den- 

sities p (half filling is given by p = 1, J = 0.35 1, t = tf). FIG - 16 - Hole densit y for a ( 192 x 2 ) *" 3 ladder at one 

quarter doping (p = 0.75) showing a pinned CDW. Here 

J' = J = 0.3* and t' = t. 
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